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We study the influence of surface tension on the shape of the conical miniscus built up by a 
magnetic fluid surrounding a current-carrying wire. Minimization of the total energy of the system 
leads to a singular second order boundary value problem for the function C,{r) describing the axially 
symmetric shape of the free surface. An appropriate transformation regularizes the problem and 
allows a straightforward numerical solution. We also study the effects a superimposed second liquid, 
a nonlinear magnetization law of the magnetic fluid, and the influence of the diameter of the wire 
on the free surface proflle. 
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I. INTRODUCTION 

The shape of the free surface of a ferrofluid in a static 
magnetic field is one of the prominent examples for a non- 
trivial interplay between magnetic and hydrodynamic de- 
grees of freedom in ferrohydrodynamics: On the one hand 
the magnetic stresses contribute to the force balance de- 
termining the surface profile whereas on the other hand 
the local magnetic field depends on this profile due to the 
magnetic boundary conditions at the surface of the mag- 
netically permeable material (1| . A particularily popular 
setup is the conical meniscus of a ferrofluid surrounding a 
vertical current-conducting wire 0, see Fig. 1. Neglect- 
ing surface tension the conical shape can be determined 
analytically from the balance between gravitational and 
magnetic force 0. 

In the present note we investigate theoretically the in- 
fluence of surface tension on the shape of the conical 
meniscus. Although surface tension is known to mod- 
ify free surface profiles in ferrohydrodynamics and its in- 
fluence has been studied, e.g., for small gaps 0,13 and 
capillaries 0, 13 to our knowledge no systematic inves- 
tigation has been done so far for the conical meniscus 
problem. Clearly, surface tension is present in all exper- 
iments in ferrohydrodynamics and a full understanding 
of its impact on a basic experiment in the field is hence 
desirable. 

As expected no analytical expression for the free sur- 
face profile can be derived when surface tension is in- 
cluded. The numerical determination can be reduced to a 
singular boundary value problem for an ordinary differen- 
tial equation which is, however, not completely straight- 
forward to solve. Still, using an appropriate transforma- 
tion, accurate solutions can be obtained. The numeri- 
cal effort is much smaller than in alternative procedures 
such as finite difference, Galcrkin, or collocation methods 
which is for example used from commercial software like 
ANSYS or COMSOL Multiphysics |1|7| 

The paper is organized as follows. In section II we col- 
lect the basic equations describing the system. Section III 
contains the analysis of a somewhat idealized situation. 
In section IV we discuss the modification of the results 



if the idealizations of section III are removed. Finally 
section V provides some conclusions. 



II. BASIC EQUATIONS 

The setup to be investigated is sketched in Fig. 1. An 
infinitely long, straight wire of radius R oriented along 
the z-axis of a cylindrical coordinate system (r, (p, z) car- 
ries an electric current /. The wire is surrounded by two 
superimposed liquids. The lower one with density qi is 
a ferrofluid with a given magnetization law M = M{H). 
The upper one has density Q2 < Qi and is a non-magnetic 
fluid. The interface between the two fluids is character- 
ized by an interface tension cri,2. For / = there is no 
magnetic field and the fiat interface between the two flu- 
ids is taken as the z — 0-plane of the coordinate system. 
For / 7^ a magnetic field 
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builds up which induces a force density in the lower 
fluid. This gives rise to an axis-symmetric conical inter- 
face parametrized by a function C(r) the determination of 
which is the central aim. Note that the magnetic fleld 
is everywhere tangential to the surface and is therefore 
independent of the detailed shape of the surface which 
makes the analysis of this case particularily transparent. 

One way to determine C,{r) is by minimizing the to- 
tal energy i^tot of the system. It is convenient then 
to use the energy of the flat interface conflguration as 
reference state and hence to minimize the energy differ- 
ence between states with a non-trivial C,(r) and C = 0. 
Equivalently one may start with the ferrohydrodynamic 
Bernoulli equation 

The total energy difference i?tot is the sum of three 
parts, the gravitational, the surface, and the magnetic 
energy. 
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FIG. 1: Sketch of the setup and definition of the main vari- 
ables. The free surface (^(r) of the ferrofluid is axis-symmetric. 



The different parts are given, respectively, by 
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Here V denotes that part of the volume occupied by the 
ferrofluid for which ^ > 0, dV is the area of the interface 
between the two fluids and jio is the permeability of free 
space. Note that in the last equality of lO we have as- 
sumed that the magnetization M and the magnetic field 
H are parallel as is the case in ferrofluids. Note also that 
H{r) is given by 

The Euler-Lagrange equation 

(6) 
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corresponding to the minimization of i^totiCl^)] is a non- 
linear ordinary differential equation for the desired inter- 
face profile C('')- For a unique solution this equation has 
to be complemented by appropriate boundary conditions. 
These are 



lim C{r) = 



and 



C'(i?) - 



1,2 
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where the prime denotes differentiation with respect to r. 
The second boundary condition results from the Young 
equation 



0-2,3 = 0-1,3 + 0-1,2 COS 61, 



(9) 



for the contact angle 6 of the fluids at the wire and the 
fact that CiR) = -tan(7r/2 - 9) (cf. Fig. 1). 



III. SIMPLIFIED MODEL 

We first analyse a somewhat simplified version of the 
boundary value problem derived in the last section. The 
simplifications are the following. We assume that the 
non-magnetic fluid is absent, g2 = 0, that the magneti- 
zation law of the ferrofluid is linear, M{H) — xH where 
X denotes the magnetic susceptibility, and that the di- 
ameter of the wire is negligible, R = 0. 

Since the magnetic field H{r) diverges for r — > so 
does the magnetic force density. In order to get a stable 
interface profile the magnetic force has to be counterbal- 
anced by gravitation and surface tension which is possible 
only if 
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This equation replaces the boundary condition JH)) in the 
present case. Consequently 0-1,3 ^^d 0-2,3 are irrelevant 
and only 01,2 remains which will be denoted simply by o 
in the present section. 

Using the simplifying assumptions we find for the mag- 
netic energy (O the expression 



E„ 



Mo X 
An 



dr 



CM 



(11) 



which allows to explicitly perform the variation in cq.JHJ. 
As a result we get the following differential equation for 
C(r-) 



M0X^1_ + C^' + rC" 1 
^^^^ 87r2 r2 (l + C'2)3/2 r 
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It is convenient to introduce dimensionless quantities by 
measuring both r and in units of 
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and a in units of giga^. Eq. H12|l then acquires the form 
1 C' + C'^ + rC" 1 



c 
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= 0. 
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This equation is easily solved for ct = yielding the well- 
known hyperbola C(r) = l/r^ P|. A perturbative solu- 
tion for CT ^ by expanding C(r) in powers of a was 
found to yield satisfactory results only for values of a 
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much smaller than those relevant in experiments. We 
therefore turned to a numerical solution. 

Eq. (|14|l together with the boundary conditions |2l and 
(|1U|I represents a singular boundary value problem of sec- 
ond kind (see, e.g., [ij). The main problem that makes 
a straightforward numerical solution inexpedient are the 
infinite intervals for r and C,. We therefore transform 
both quantities according to (cf. also Fig. 1) 
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and describe the interface profile by p = pi'fi)- The 
boundary conditions Q and H1U|I then transform into 



respectively which are much more convenient for the sub- 
sequent numerical solution. Substituting the derivatives 
in ifTHl according to 
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(16) we find as differential equation for p{Lp) 
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Rewriting this equation in the form p'^f) = 
f{(p,p{(p),p'{(p)) the boundary value problem l(T^ . l(TB|) 
can now easily be solved numerically using, e. g., the 
nonlinear finite-difference method described in [12| . As 
initial guess for the solution which is needed in this 
procedure the transformation of the analytical solution 
Co{r) = l/r^ for cr = given by po{(p) = cos{ip) ^ ta,n{ip) 
may be used. 

Fig 13 shows the resulting interface profile for param- 
eter values corresponding to the ferrofluid EMG909 of 
Ferrotec |0 . For comparison the shape for cr = is also 
shown. It is clearly seen that the free surface profile is 
markedly modified by the influence of the surface ten- 
sion. As expected the inclusion of surface tension makes 
the profile narrower since an additional force pointing 
radially inward builds up. 
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IV. MORE REALISTIC SITUATIONS 

The above solution of the idealized problem forms a 
convenient starting point for the discussion of the in- 
fluence of those feature in the original setup that were 
neglected in the previous section. 

A simple modification is the case in which the ferrofluid 
is superimposed by a non-magnetic liquid. Then the den- 
sity gi has to substituted by the density difference gi — Q2 
and the surface tension a is to be replaced by the interface 
tension ai^2- The gravitational contribution to the to- 
tal energy gets reduced and consequently larger displace- 
ments from the flat interface are to be expected. The 
situation can be analyzed quantitatively without further 
effort by mapping it to the case analyzed in the previous 
section. In fact the interface proflle is again determined 



FIG. 2: (a) Free surface profile C,{r) as determined numerically 
for i? = and / = 30A for the ferrofluid EMG909 with 
parameters gi = 1120 kgm~^, a = 0.0259Nm"\ and x = 0.8 
(full line). Also shown is the result for a — and otherwise 
identical parameters (dashed line) . (b) Same profiles as in (a) 
in terms of the transformed function p{(p). 



by eq. H14|l with only the dimensionless units and the 
value of a being modifled. As an explicit example we 
show in Fig.|3Ia) how the interface proflle of Fig. [21 gets 
modifled for g2 = lO'^kgm^^ and tTi_2 = cr. 

Deviations form the linear magnetization law which 
become relevant in particular for small values of r where 
the magnetic field gets strong can also be dealt with. 
An improved approximation is provided by the Langevin 



4 



relation 0, 



M{H) = M, coth(a 



a ■ 



mH 
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with the saturation magnetization of the fluid Mg, the 
magnetic moment m of a single ferromagnetic parti- 
cle, the Boltzmann-constant fcs and the temperature 
T. The zero-field susceptibility is then given by x = 
mM^ / {ik-QT) . The density of the magnetic energy can 
again be determined analytically 



dH ■ M(H) = ^ ^ \n . (21) 



This expression gives rise to a somewhat more compli- 
cated equation for C('")i however, the overall procedure 
of section III remains valid. As shown in Fig. Ota) the 
modification of the result for linear magnetization law is 
usually small in accordance with the fact that the frac- 
tion of the ferrofluid volume for which the field is large 
is rather small. 

Finally, to elucidate the influence of a non-zero radius 
of the wire we consider the problem for R > but with 
g2 — and a linear magnetization law. The main qual- 
itative difference is that instead of Hl()|l we have to use 
as boundary condition for small r. This implies that 
C(r) is not singular at the lower boundary. Hence al- 
ready the simple transformation r = 1/m is sufficient to 
get a numerically tractable boundary value problem for 
the function C(u). As initial guess a linear dependence 
Co('") = Am may be used. 



Note that A fixes the value of C{r) at r = R whereas 
the boundary condition © involves the derivative C'{r) 
at r = i?. Therefore A or equivalently C{R) has to be 
modified until the required value for C{R) is obtained. 
This procedure leads to a monotonic relation between 
the prescribed contact angle 9 and C{R), cf. Fig. I^l^c). 
With the help this value for C(^) the complete surface 
profile C(^) can then be determined. Fig. Otb) shows 
a collection of surface profiles with different C(^) ^'^d 
corresoponding contact angles 6. In the limiting case 
R ^ 0, the surface profile is independent of the contact 
angle 9 as demonstrated in Fig. Eld) for a sequence of 
profiles with different values of R and a fixed contact 
angle 9 = 45°. 



V. CONCLUSIONS 

In the present investigation we have determined the 
free surface profile of a magnetic fluid surrounding a ver- 
tical current carrying wire with special emphasis on the 
effects of interface and surface tension respectively. We 
have found that for experimentally relevant parameter 
values there is a strong influence of the surface tension 
giving rise to a more slender profile as compared to the 
well-studied case without surface tension P| . A more gen- 
eral magnetization law including saturation as well as the 
prescription of the contact angle at the wire changes the 
profile only in the vicinity of the wire. A superimposed 
non-magnetic liquid changes the overall scale of the pro- 
file due to the reduction of the gravitational energy. 
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FIG. 3: Modifications of tfie idealized results of section 3. If not explicitly stated otherwise the parameters are the same as in 
Fig. |5| The dashed line shows always the same interface profile as the full line in Fig. |5| (a) The dotted line shows the interface 
profile if a non-magnetic fluid with density g2 — 10"^ kgm^"^ is superimposed to the ferrofiuid. Also in (a), the solid line shows 
the interface profile for a ferrofiuid with magnetization law of Langevin type with Ms = 16kAm~^ and the same zero-field 
susceptibility x ^ in Fig- ID (b) Dependence of the interface profile on the contact angle 9 for non-zero radius R = 0.5 mm. 
(c) Relation between the contact angle 9 and the height of the interface at the wire C,{R). Marked points correspond to the 
profiles shown in (b). (d) Sequence of surface profiles with fixed contact angle 9 = 45° for decreasing diameter of the wire, 
7? = 2; 1.5; 1; 0.5; 0.2 mm. For _R — » the profile of Fig. |21is almost reproduced. 



